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Abstract 

We propose coupled evolution equations for the thickness of a liquid film and the density of an adsorbate 
layer on a partially wetting solid substrate. Therein, running droplets are studied assuming a chemical 
reaction underneath the droplets that induces a wettability gradient on the substrate and provides the driving 
force for droplet motion. Two different regimes for moving droplets - reaction-limited and saturated regime 
- are described. They correspond to increasing and decreasing velocities with increasing reaction rates 
and droplet sizes, respectively. The existence of the two regimes offers a natural explanation of prior 
experimental observations. 
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It is well known that a drop of liquid can start to move if it is exposed to an external gradient. 
For instance, a drop of liquid immersed in another liquid that is placed in a temperature gradient 
will move towards the higher temperature due to Marangoni forces caused by surface tension 
gradients [1]. Alternatively, a temperature gradient [2] or a wettability gradient induced by a 
chemical grading of the substrate [3, 4] cause drop motion. 

Drops may also move in initially homogeneous surroundings. This is possible because such 
active drops change their surrounding and produce an internal gradient that drives their motion. 
One example are drops immersed in another liquid that contain a soluble surfactant undergoing an 
isothermal chemical reaction at the surface of the drop [5, 6]. Recent experiments found a differ- 
ent sort of chemically driven running droplets on solid substrates. There, the driving wettability 
gradient is produced by a chemical reaction at the substrate underneath the drop [7-9]. 

In these experiments, a small droplet of solvent is put on a partially wettable substrate. A 
chemical dissolved in the droplet starts to react with the substrate resulting in the deposition of a 
less wettable coating. The substrate below becomes less wettable than the substrate outside the 
droplet. Eventually, the symmetry is broken by fluctuations and the drop starts to move away from 
the 'bad spot', thereby changing the substrate and leaving a less wettable trail behind (see Fig. 1). 
Similar phenomena can be seen in reactive (de) wetting [10], in camphor boats [11], and in the 
migration of reactive islands in alloying [12]. 




substrate 

FIG. 1 : Sketch of a right moving droplet driven by a chemical gradient. 

From a simple theoretical argument [13], an implicit equation for the velocity v of the droplet 
was derived [13]: v = C tan6'*(l — exp(— rL/f )), where r is the reaction rate, L the droplet's 
length and C a constant. The dynamic contact angle 6* is assumed to be same at the advancing 
and receding ends of the droplet, i.e. the droplet profile is approximated by a spherical cap. 9* is 
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given by cos^* = (cos^" + cos 91) /2, where the static contact angles at the advancing edge 9^ 
and at the receding edge 91 > 9^ are different due to the chemical gradient. The expression for the 
velocity is found assuming a first order reaction on the substrate yielding chemical concentrations 

= and a,. = 1 — exp{—rL/v) at the respective ends of the droplet. The expression for the 
velocity predicts a monotone increase of the droplet velocity with the droplet length L and the 
reaction rate r, in line with early experimental observations [8]. Recent experiments [9] show the 
opposite trend; the velocity decreases with increasing drop sizes and reaction rate. 

In this Letter we propose and analyze a dynamical model for self-propelled running droplets. 
It consists of coupled evolution equations describing the interdependent spatiotemporal dynamics 
of the film thickness h and a chemical concentration a that changes the substrate wettability. This 
model is capable of reproducing both experimentally found regimes. In particular, we identify two 
distinct regimes of running drops separated by a maximum of the velocity. For small droplet size 
or reaction rate, the chemical gradient in the drop is limited by the progress of reaction. In contrast 
for large droplets and fast reaction, the chemical concentration at the receding end saturates at a 
maximum value. Below, we show that the velocity of reaction-limited droplets increases, while 
the velocity of saturated droplets decreases with increasing reaction rate and drop size. In addition, 
we find that the dynamic contact angles at the advancing and receding edges of the droplets differ 
substantially. If we take into account fast diffusion of the chemical substance, the droplet velocity 
in the saturated regime decreases much faster and only stationary drops are found at large reaction 
rates. The velocity of the droplets may increase [8] or decrease [9] with their volume. 

A liquid layer on a smooth solid substrate is modelled by an evolution equation for the film 
thickness profile h{x, t) derived from the Navier-Stokes equations using long-wave approximation 
[14] 



dfh = -V 



^V{^Ah + Il{h)) 
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with 7 and -q being the surface tension and viscosity of the liquid, respectively. The disjoining 
pressure n(/i) accounts for the wetting properties of the substrate [15]. We use the form Il{h) = 
2Sadl/ + (Sp/l) exp[((io — h)/l] [16, 17], where Sa and Sp are the apolar and polar components 
of the total spreading coefficient S = Sa + Sp, respectively, = 0.158 nm is the Born repulsion 
length and I is a correlation length [16]. We choose Sa > and Sp < thereby combining a 
stabilizing long-range van der Waals and a destabilizing short-range polar interaction. The latter 
contains the influence of the coating and crucially influences the static contact angle [16]. Such a 
model allows for solutions representing static droplets with a finite mesoscopic equilibrium contact 



angle sitting on an ultrathin precursor film. 

To account for the varying wettability caused by the chemical reaction we let the polar part 
of the spreading coefficient Sp depend linearly on the coating density, i.e. on the concentration 

a{x, t) of a chemical species adsorbed at the substrate: 

5p = ^°(l + ^)<0. (2) 

The equilibrium contact angle 9e is given by cos ^e = '5'/7 + l[16]. This implies that 9e increases 
with a, i.e. the coated substrate is less wettable. The constant Sp / g' defines the magnitude of 
the wettability gradient. Note, that Eq. (2) corresponds to the linear relation between cos 9e and a 
assumed in [8, 9, 13]. 

The evolution of the chemical concentration on the substrate is modelled by a reaction-diffusion 
equation for a{x) 

dta ^ R{h,a) + d' Aa, (3) 

where the function R(h, a) describes the reaction that changes the wettability of the substrate and 
the second term allows for a (usually small) diffusion of the chemical species along the substrate. 
The main results, however, are obtained without diffusion. As reaction term we choose 

R{h, a) = r'e{h - ho) (l - . (4) 

The time scale of the reaction is defined by the rate constant r'. It is assumed that the reaction at 
the substrate occurs only underneath the droplet; this is modelled by the step function Q{h — ho). 
The value of ho is chosen slightly larger than the thickness of the precursor film. The chemical 
concentration of the coating saturates at a value ctmax, because the reaction is assumed to be fast 
enough to produce a complete adsorption layer. Changes in details of the reaction term like the 
replacement of the step function by an explicit proportionality to h do not affect our results qual- 
itatively. The physical motivation behind the actual choice is a fast equilibration of the reactant 
concentration within the moving droplet, caused by diffusion and convective motion within the 
droplet. 

We rewrite Eqs. (1) to (4) by introducing scales 3jr]/lK^, ^Ij/k, a^ax and / for t, x, a and h, 
respectively. Then, k — {\Sp\/l) exp(do/0 is the dimensionless spreading coefficient. Defining the 
dimensionless reaction rate r = Sr'^rj/lK^, diffusion constant d — Sd'Tj/nl"^, wettability gradient 
g — g'/aiaax and ratio of the effective molecular interactions b — 2Sadl/l^K,, we obtain from 
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FIG. 2: Profiles of the droplet h{x, t) ((a) and (b)) and concentration a of coating ((c) and (d)) in the two 
qualitatively different running droplet regimes, (a) and (c) Reaction-hmited regime at r = 0.0001 moving 
with velocity v = 0.026. (b) and (d) Saturated regime at r = 0.001 with v = 0.032. The streamlines plotted 
in (a) and (b) indicate the convective motion inside the droplet in the comoving frame [19]. The remaining 
parameters are = 1.0, d = 0.0, h = 0.5 and the droplet volume is 30000. 

Eqs. (1) to (4) the dimensionless coupled evolution equations for the thickness profile h and the 
substrate coverage a 



dth = -V ■ yi'v 

dta = rQ{h — ho) {1 — a) + d Aa 



a 
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In the one-dimensional version of Eqs. (5,6), we calculate by use of continuation techniques [18] 
two-dimensional running droplets moving with constant speed. This is achieved by switching to 
the comoving frame x — vt and imposing appropriate boundary conditions [19]. 

Fig. 2 shows the constant profiles of such moving droplets ((a) and (b)) together with the corre- 
sponding concentration profiles for a ((c) and (d)) for two qualitatively different regimes (promi- 
nently visible in the a profiles). In ((a) and (b)) the streamlines in the comoving frame indicate the 
convective motion in the droplets. 



In Figs. 2 (a,c) the limiting factor for the motion is the small reaction rate. The reaction time 
is long compared to the viscous time scale, i.e. the droplet passes too fast for the coverage a to 
approach its saturation value. In the case = 0, we can estimate the value of a at the receding end 
to be = 1 — exp(— rL/i;) where L is the length of the droplet. We call this regime reaction- 
limited. In Figs. 2(b,d) the reaction at the receding edge has reached its saturating value, i.e. 
Or ~ 1 and rhjv >> 1. After the droplet has passed, the concentration a stays at the saturation 
level. We call this regime saturated. One can not exactly predict where the transition between the 
regimes occurs, because the velocity of the droplet and the driving wettability gradient determine 
each other dynamically. We found that in both regimes the advancing contact angle is substantially 
smaller than the receding one. This contradicts the assumption of identical dynamic contact angles 
in the approximate theory [13]. 
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FIG. 3: Moving drops in three dimensions. Shown are the film thickness profiles in the (a) reaction-limited 
and (b) saturated regime, (c) and (d) represent the corresponding contour lines. The velocity converged to 
V ^ 0.018 and v ^ 0.044, respectively. The direction of movement is indicated by arrows, the parameters 
correspond to Figs. 2 (a,c), respectively, and the droplet volume is 7.5 x 10^. Only a part of the computational 
domain size of 500 x 500 is shown. 
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FIG. 4: Characterization of running droplets stationary in a comoved frame. Shown are (a) the velocity v in 
dependence of the reaction rate r (logarithmic scale) without (d = 0, solid Une) and with (d = 1.0, dotted 
line) diffusion.; (b) the tangents of the advancing 9a (solid line) and receding 6^ (dashed Une) dynamic 
contact angles in the reaction-limited regime; (c) the absolut value of the difference of tan 9a (soUd line) 

and tan 6r (dashed line) and the tangents of the static contact angle 9e{a) at the corresponding a values. 
Also in (c) only the reaction-limited regime r < is shown where always 9a > 9e{aa) and 9r < 9e{oir). 
The remaining parameters are as in Fig. 2. 

We have tested that the above model also allows for stationary moving droplets in three di- 
mensions (Fig. 3) by integrating Eqs. (5) and (6) in time starting from adequate initial conditions 
[20]. After a transient the droplets approach constant shape and velocity. We show resulting 
droplets for the reaction-limited (Figs. 3 (a,c)) and the saturated (Figs. 3 (b,d)) regime. Both drops 

have oval shape with an asymmetry between advancing and receding part that is stronger in the 
saturation-limited regime due to the larger wettability gradient. 

A closer analysis of the one-dimensional version of Eqs. (5,6) gives insight into the properties 
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of the running droplets and their existence and stability conditions. Fig. 4(a) shows the droplet 
velocity in dependence of the reaction rate r without and with diffusion. The dependence of the 
velocity on reaction rate r for the running droplets is non-mono tonic. There exists a value 
where the velocity is maximal. The case r < [r > r^] defines the reaction-limited [saturated] 
regime. 

The velocity increase for r < is easily explained by the increasing value of (Xr at the back 
of the drop caused by the faster reaction. The decreasing velocity for r > originates from 
a slow increase of ol in the contact region at the front of the drop, while a at the back remains 
constant at its saturation value {ol,. ~ 1). We have also studied the case where diffusion of a. 
plays is sizable. There, the rapidly produced a diffuses to the substrate in front of the droplet 
and diminishes the driving gradient between front and back of the drop. The velocity curve is 
similar for small r but decreases even faster beyond the maximum and drops to zero at large 
enough r; there, no running droplets are found anymore. The latter phenomena is reminiscent 
to observations of moving spots in reaction-diffusion systems of activator-inhibitor-type [21-23]. 
Therein, a transition from moving to stationary droplets (drift pitchfork bifurcation) occurs when 
the inhibitor is produced to fast. 
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FIG. 5: The transition between (i) reaction-limited and (ii) saturation-limited regime is most clearly seen in 
the dependence of velocity v on drop size, i.e. droplet volume V . Curves are plotted for different reaction 
rates r (see legend). The remaining parameters are as in Fig. 2. 

Figs. 4(b) and (c) give, respectively, the tangents of the measured advancing and receding dy- 
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namic contact angles (tan^a/^) and their differences Atan^a/^ — tan 9a/r — tan^e(«a/r-) from 
the static contact angles as a function of the droplet's velocity in the case d — 0. All contact 
angles given correspond to the mesoscopic contact angles measured at the inflection points of 
the computed drop profiles. To determine the static contact angle we take the value of a at the 
inflection point to calculate a space independent polar spreading coefficient in Eq. (5). Fig. 4 (b) 
shows that the two dynamic contact angles may differ quite strongly. However, a comparison of 
Figs. 4 (b) and (c) shows that the difference between the static and the dynamic contact angles is 
much smaller. 

Experimentally, an important measure is the dependence of droplet velocity on droplet length 
or volume [8, 9]. We show in Fig. 5 the dependence of v on the droplet volume. For small 
reaction rates r we find that the velocity increases with drop length. The increase becomes more 
pronounced for larger r. This corresponds to the reaction-limited regime and the experiments in 
[8]. Upon further increase of r, the curve shows a velocity maximum at an intermediate drop size. 
For bigger drops, the velocity decreases monotonic similar to the observations of [9]. Inspection 
of the corresponding profiles suggests that the regime of decreasing velocity is again the saturated 
regime. 

To conclude, we have modelled chemically driven moving droplets by a dynamical model con- 
sisting of coupled evolution equations for the film thickness profile and the concentration of a 
surface coating that is produced underneath the droplet. This model exhibits moving droplets with 
constant speed and shape in two- and three dimensions. A systematic analysis of the velocities 
of two-dimensional droplets allows us to distinguish the reaction limited and saturated regimes, 
which are characterized by an increase resp. decrease of droplet velocities with increasing droplet 
size or reaction rate. The behavior in the saturated regime indicates that the velocity of the droplets 
in our model depend not only on the values of the chemical coating at the advancing and receding 
ends of the droplet, but are in addition determined by the concentration profiles of a near the con- 
tact lines. The two regimes offer a natural explanation of the conflicting experimental observation 
in [8, 9]. 
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